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ABSTRACT 

We present a simple physical model for populating dark matter (DM) halos with Lyman Alpha 
Emiiters (LAEs) and predict the physical properties of LAEs at z ~ 3 — 7. The central tenet of this 
model is that the Lya luminosity is proportional to the star formation rate (SFR) which is directly 
related to the halo mass accretion rate. The only free parameter in our model is then the star- 
formation efhciency (SFE). An efficiency of 2.5% provides the best-fit to the Lya luminosity function 
(LF) at redshift z = 3.1, and we use this SFE to construct Lya LFs at other redshifts. Our model 
reproduce the Lya LFs, stellar ages, SFR« 1 — 10 Mq yr~^, stellar masses ~ 10^ — 10^ Mq and the 
clustering properties of LAEs at z « 3 — 7. We find the spatial correlation lengths 3 — 6 ft.^^ Mpc, 
in agreement with the observations. Finally, we estimate the field-to-field variation « 30% for current 
volume and flux limited surveys, again consistent with observations. Our results suggest that the star 
formation, and hence Lya emission in LAEs is powered by the accretion of new material, and that 
the physical properties of LAEs do not evolve significantly over a wide range of redshifts. Relating 
the accreted mass, rather than the total mass, to the Lya luminosity of LAEs naturally gives rise to 
the duty cycle of LAEs. 

Subject headings: galaxies: high-redshift — galaxies: Lyman alpha emitters — galaxy: luminosity 
function — galaxy: clustering — galaxy: correlation length 



1. INTRODUCTION 

Lyman-a emitting galaxies (LAEs), which seem to 
be powered by star formation activity, are selected 
on the basis of strong Lya emission li ne, irrespec- 
tive of other galaxy pr o pertie s (e.g. | Cowie 



2001 



tive o r other galaxy pr o pertie s (e.g. 
1991 iRho ads et all l2000, (20p; iFvnbo e t al . . _ , . 
Aiiki et all 2003; Matsuda et al. 2005; Tani guchi et al l 
20051 iGa wiser et al. 2006; Shimasaku et all \20oS. 



Tapken et al. 2006.: Mur avama e t al. 2007 ; Nilsson et al.l 
2007t lOuchi et all 12008^ However, the high-redshift 



galaxies selected on the basis of this one property are 
reasonably uniform in some of their other properties^. 
For example, the inferred stellar masses of LAEs at z 
< 5 a r e typically sma ll , ~ 1 0^ — lO^Mp (IGawiser et al.l 



2006t iPirzkal et al.l 120071: iFinkelstein et al 
Pentericci et al. I I2009D . and often have 



[20071 : 

large Lya 

equivalent width ind i cating a young stellar population 
(jMalhotra fc Rhoadi l2002l ) . which is also sup ported 
jthe blue colo r of these galaxies ([Venemans et al.l 



2005; Finkclstein e t al.ll2007l 120081 iGronwall et"all T 
Pirzkal et al.. .2007) . especially at high redshifts. AGNs 



are ruled out as sources of strong Lya er nission in LAEs 
due to non-detection of X-ray emission (Mal hotra et al.l 
[2003; Wang et al. 2004; Gawiser ct al. 2007) and lack of 
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high i onizat ion lines in the optical spectra ([Dawson et al.l 
|2C)04 [2007D . LAEs hav e a in o derate SFR »5 - 8 Mq, 
yr-i (e.g. iPirzkal et all [20071 : [Taniguchi et~al[ [2005f) 
and spatial correlation lengths of » 3-5 ..-M^PC; 
albeit with a substanti al uncertaintv fOuchi et al.l [20031 : 



[Kovac et al. 2007; Ga wiser et al.|[200 

Despite the increasing number of LAE observations, 
theoretical understanding of LAEs is still in early 
stages primarily due to poor understanding of phys- 
ical properties including star formation rate (SFR), 
stellar initial mass function, Lya escape fraction, and 
the duty cycle of Lya phase of star-forming galax- 
ies at high redshifts. There have been several the- 
oretical studies of LAEs based on cosmologica l sim- 

"2001 

2007; Nagaminc et all 
(e.g. iLe Delliou et al.l 



iTasitsiom il [^6: Shimizu et al. 
'2008') , "semi-analytical mo dels 



Haiman & Soaans 
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L999I: [Diikstra et al [[20071 [Mao et all 


20071: [Stark et al. 




20071 [Fernandez & Komatsul [2008D 



that relate the total halo mass to the Lya luminosity of 
LAEs. Such a linear relationship between the halo mass 
and the Lya luminosity often leads to an overprediction 
of the number density of LAEs. To reconcile the mass 
distribution of halos to the LF one needs to either assume 
a small escape fraction of Lya photons (which would fail 
to account for the la rge Lya EWs observed) or intro- 
duce a duty cycle (e.g. [Stark et al.|[2007l : [Nagamine et al.l 
[2008| ) which adds another parameter to the models. In 
addition, the complexity and large number of variable pa- 
rameters in many models motivates the development of 
a simple model that can reproduce the LAE observations 
over a wide range of redshifts. Such a model would be 
particularly useful in understanding the nature of LAEs 
observed at high redshifts. 

In this paper we present a physical model to populate 
DM halos with LAEs in a cosmological simulation, and 
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predict the abundances & physical properties of LAEs 
at z « 3 — 7. This model differs from many of the ear- 
lier studies in a major way, namely in that we relate 
mass accreted, as opposed to the total halo mass, to the 
Lya luminosity. Mass accretion onto halos via smooth 
infall and accretion due to mergers of a specific mass 
ratios has been s hown to have distinctly differe nt clus- 
tering behaviour (jScannapieco fc Thackeil [2003D . How- 
ever, in this model we do not distinguish between the 
smooth accretion and the accretion due to mergers. In 
other words, in our model the LAEs are undergoing an 
episode of star-formation which is driven by accretion of 
fresh material onto the halos, independent of whether the 
accretion is due to mergers or if it is via smooth infall. 
While there is no direct observational evidence showing a 
relation between the bary ons accreted and the Lya lumi- 
nosit y, recent studies (e.g. iDekel et"aLll2009l : iKeres et"al] 
l2009f) have shown that such cold accretion of new mate- 
rial can drive star formation in galaxies, especially those 
hosted by less massive halos. In this model we assume 
that the LAEs do not contain large amounts of dust, 
and hence most of the hydrogen ionizing photons will 
be absorbed, while most of the Lya photons can escape 
(jGawiser et al.l 120061: iKobavashi et all l2007f) . This as- 
sumption yields large Lya equivalent widths even with- 
out appealing to metal-free Population III stars, whose 
contributions ar e constrained by non-det e ction of the 
Hell (1640) line (|Dawson et al.ll2004 [20071 iNagao et all 
120081 ). In such condit ions, the Lya lin e becomes a direct 
measure of the SFR (! Kennicutti[l98l . which is propor- 
tional to the accretion of fresh gas onto the galaxy. The 
constant of proportionality (the star formation efficiency) 
between accretion rates and Lya luminosity is the only 
free parameter in our models. 

This study is organized as follows: In §2 we give a 
detailed description of our physical model. We describe 
the DM simulation parameters, and how we generate DM 
halo catalogs in §3. In §4 we first construct Lya luminos- 
ity function(LF) using model LAEs at z « 3 and compare 
it to the observations to find the best-fit model parame- 
ter, and then use this best-fit parameter to construct Lya 
LFs at z > 3. In §5 we derive the physical properties of 
LAEs using our best-fit model and then investigate the 
evolution of Lya LF. We study the large scale structure 
of model LAEs, and study the redshift evolution of corre- 
lation lengths of LAEs in §6. We summarize and present 
conclusions in §7. 

2. PHYSICAL MODEL FOR LYMAN ALPHA EMITTERS 

Our model is motivated by the idea that Lya 
emission in LAEs is ass ociated with star-formation 
(jPartridge &: Peeble£lll967| ) from rapid accretion of new 
material on to the DM halos. This new material pro- 
vides fresh fuel to the system driving the star formation 
(jKeres et al.ll2009( ) a nd consequently increasing the L ya 
line flux in galaxies (jKennicuttlliOSa iHu et alJIloM l. 

In our model we populate each DM halo with an LAE, 
and assig n to it Lyo: luminosity (Llvo) p roportional to 
the SFR fKen nicuttllTQSl iHu et al.lll999l ) using the fol- 
lowing equation; 

SFR 



L 



Lya 



1 X 10*2 



.,r _i ergs . (1) 

Here L^ya is the intrinsic luminosity of an LAE. The ob- 
served Lya flux will also depend on the escape fraction of 



the Lya photons (/i^c")- Moreover, equation (1) implic- 
itly assumes an escape fraction near zero for the ionizing 
continuum photons, whose absorption is required to pro- 
duce the Lya emission line. 

While the escape fraction of ionizing Lyman contin- 
uum photons (A < 912 A) is not very precisely known, 
several studies have shown that an escape fraction of 
only a few percent is sufficient to meet the observ ational 
constraints on the reionization epoch (e .g. Wood fc Loeb 
12000: Hansen fc Oh 2006|; iRazoumov fc Sommer-Larsen 
|2006 : _Gncdin_ct_al. 200^. In addition, observations are 
also generally consistent with small esc ape fractions of 
Lyma n continuum photons both loc allv ([Leitherer et al.l 



19951) and at high redshifts (e.g. ISteidel et al.l |2001| 
Hans en fc OU 120061 : IShaplev et al.l l2006t llwata et all 
,20091 

The escape fraction of Lya photons {feJi") is ob- 
served to be large wit h f^^f' ft! 1 (jGawiser et al.l 120061 : 
IKobavashi et al.ll2007D likely causing the observed large 
Lyg equivalen t widt hs fEWs) (iMalhotra fc Rhoadsll2002l : 
iHansen fc OhI l2006l : iFinkelstein et al.ll2007l] . However, 
the semi-analytic model of lLe DeUiou et ail (j2006t ) pre- 
dicts a much smaller value of /i^^" = 2%. For simplicity, 
in our model we approximate f^Jl"' and Lyman contin- 
uum photons as unity and zero, respectively. Thus all 
the Lya photons produced in LAEs escape to be ob- 
served while none of the ionizing photons escape from 
the galaxy. Small deviations from these assumptions will 
not affect our results significantly. 

As stated earlier, in our model we assume that the 
rapid accretion of new material on to the DM halos 
causes star formation in LAEs. We estimate the SFR 
(mass in accreted gas, AMgas converted to stars in 
unit time) in LAEs by converting baryonic mass ac- 
creted {AM},) by DM halos, adopting a constant ratio 
of baryons to the dark matter, over a short time scale, 
tLya- This time scale (« 30 Myrs ) is broadly simi- 
lar to the stellar population ages of most Lya galaxies 
(e.g. [Pirzkal et al. 2007; Finkclstcin ct al. 2007, 2008), 
the lifetimes of OB associations, and the dynamical time 
expected for L ya galaxies based on their measured sizes 
(jPirzkal et al.l f2006. 20 Q1). A sim i lar tim e scale (« 70 
Myrs) was used by iShimizu et al.l ([20071) to match the 
morphology of large scale structure of LAEs by varying 
the amplitude of density fluctuations in galaxy formation 
models. 

Thus, 



SFR^Ux 



{AM,as , ^^^^ 



\ ^Lya 



\ ^Lya 



= hxAh . (2) 



where /* is the star formation efliciency. In the above 
equation we have assumed that AMgas is same as the 
baryonic mass accreted (AAf^) by the DM halos, and 
since our simulation contained only DM particles, we 
use the universal ratio of baryonic and DM densities i.e. 
AMb = (fib/riDivi) X AMdm, where fifc and fioM are the 
baryonic and DM density parameters and AMjjm is the 
dark matter mass accreted by the DM halos. 

Finally, the total mass in young stars in an LAE, is 
estimated using 



~ SFRxtLya = h^MbXtLy 



/*x-|^xAA/dm 



DM 



(3) 
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Fig. 1. — Evolution of Lya LFs at redshifts z 3 — 7. The dotted lines show results from our model and the symbols with error bars 
are the observational data, (a) The best-fit model Lya LP at z = 3.1 yields a star formation efficiency of 2. 5%. We use this star f ormation 
efficiency to constru ct model Lyg LFs at z=4.5 , 5.7, and 6.6 (b-d). T he references for the data: z = 3.1 IIGronwall et al.ll2007l) . z = 4.5 
dPawsori et al.|[2007^ . z = 5.7 l lOuchi et al.|[200l) . and z = 6.6 JKashik awa et al. 2006). 



This corresponds to the mass of stars younger than 30 
Myrs which contribute to the Lya and UV continuum 
emission, which is more easily measured than the total 
stellar mass. The only unknown variable in all of the 
above equations is f^,, the only free parameter in our 
model. 

3. SIMULATION & HALO CATALOGS 

We generate the dark matter (DM) halo catalogs using 
an N-body DM cosmological simulation code GADGET2 
(|Springe]| l2005| ). We generated the initial conditions 
for the simul ation using 2nd-ord e r Lagrangian Perturba- 
tion Theory (jCrocce et al.ll200^ iThacker fc Couchman I 
120061) . In this simulation we use 1024^ DM particles in a 
comoving volume of (102 Mpc)'^, a volume greater than 
a typical LAE survey. Each DM particle has a mass 
« 2.7x IO'^A/q/i"^ Using Friends-of-Friends (EOF) halo 
finder, we identify DM halos that contain 100 or more 
DM particles. This corresponds to a minimum halo mass 
« 2.7 X lO^MQh^^. We then generate catalogs, for red- 
shifts from z = 10 to z = 3, which contain positions of 
halos, their DM mass, and unique IDs of each individual 
particle that belongs to a given halo. These unique parti- 
cle IDs are later used to track halos between two epochs. 
Throughout this work we assumed a flat ACDM cosmol- 
ogy with parameters f^,„=0.233, f7A=0.721, f^b=0.0462, 
ft,=0.71, CT8=0.817 where flm, ^A, ^b, h, and cor- 



respond, respectively, to the matter density, dark en- 
ergy density, and baryonic density in units of the crit- 
ical density, the Hubble parameter in units of 100 km 
s^^ Mpc""'^, and the RMS density fluctuations on the 8 
Mpc scale, in agr eement with WMAP-5 year results 
(|Hinshaw et al.ll2009[ ). 

4. LYMAN ALPHA LUMINOSITY FUNCTION 

We now construct the Lya luminosity function (Lya 
LF), the number of LAEs per unit volume in a given 
luminosity bin. First, we calculate the total dark matter 
mass accreted {AM^m) by each DM halo at z — 3.1 
during an interval « 30 Myrs (equals to t^y^ in equation 
2). To calculate AM^m we track each halo, using unique 
ID associated with particles in a given halo, between two 
epochs separated in time tLy^- In general, we expect 
every halo to accrete more mass with time. However, due 
to group finding noise, we find that some halos lose mass 
with time (negative mass accretion). In other words, the 
mass accretion by some halos is not real but results from 
the simulation noise. The main reason for this noise is 
the way halos are identified in any DM simulation. In the 
following paragraph, we discuss how we take into account 
this noise resulting from the halo finding technique. 

In our DM simulation, we use a EOF halo finder that 
identifies DM halos. In FOE, all the particles within a 
certain linking length from each other are linked to a 
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Fig. 2. — Accreted mass function and halo mass functions. 
Solid lines show accreted mass functions at ^ = 3.1 (violet), 
z = 4.5(blue), z = 5.7 (green) and z = 6.6(red). The dashed 
blue line shows the dark matter halo mass functions at ^ = 4.5 
to compare with the corresponding accreted mass function. 
The vertical dotted lines enclose the region of observed Lyo 
luminosities in LAEs. 

halo, independent of whether a given particle is gravita- 
tionally bound to a halo. Thus, associating a particle to a 
halo based on the linking length give rise to some uncer- 
tainty in halo mass (in this case AM^m)- To determine 
how many halos have real accretion (rather than spurious 
apparent accretion due to uncertainty in particle associ- 
ation to a halo by the halo finder), we first construct a 
histogram of AM^im including the halos with negative 
AMjjM- We then subtract the halo counts in negative 
AMuM bins from the corresponding counts in the posi- 
tive AMdm bins. This procedure compensates for halos 
that show accretion just due to random nature of FOF 
halo finder. The residual halos with positive accretion 
rates are then considered for constructing Lya LFs. 

Next, we convert the accreted mass bins to the Lya lu- 
minosity bins using equation 1 to yield Lya LF. We then 
compare this luminosity function with the observations 
at z = 3.1 ([Gronwall et al. 2007), and get the best-fit 
model by varying the star formation efficiency (/^) to 
yield the least reduced chi-square (chi-square per degree 
of freedom) given by 



N 



-1^ a. 



odel 



2 

model 



(4) 



'obs 



where iV, Nmodei & Nobs are the number of observed data 
points, LAE counts from model, and the observed LAE 
counts in each bin, respectively, and the Poisson errors 
are given by amodei = V Nmodei and aobs = VNobs- Fig- 
ure l(top left) shows the best-fit model Lya LF (dotted 
line) at z = 3.1 . The symbols are the observations from 
iGronwall et al.l ()2007f ) shown with 1-sigma error bars. 

Lastly, we use the best- fit model parameter i.e. the 
star formation efficiency at z = 3.1 to construct the 
model Lya LFs at z > 3, and then compare these LFs 
with the observations at z =4.5, 5.7, and 6.6. Fig. 1 
shows the differential Lya LFs from our model (dotted 
lines) and observation s (filled circles) at redshifts z = 3.1 
(IGronwall et al.ll2007D. z = 4.5 (Dawson et al.ll2007D. z = 
5.7 (jOuchi et al.ll2008[ ). and z = 6.6 (jKashikawa et all 
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Fig. 3. — The DM accretion rate of halos as a function of halo 
mass at z fa 3 — 7. The solid line is the least-square fit to the 
median mass in each AM/dt bin. The slope of the lines is nearly 
constant fx 0.8 — 0.9 over all redshifts. 



l2006f ). We have rebinned the observational data for 
z = 3.1 and z=6.6 data so as to make the bin size uni- 
form at all redshifts. Wiggles seen, especially in z = 4.5 
model Lya LF (Fig. 1 top right) are probably due to 
simulation noise. The best-fit model for z = 3.1 Lya LF 
yields a star formation efficiency of 2.5%. Corresponding 
to this star formation efficiency, the values between 
our model and the observed Lya LFs are 0.5, 0.8, 1.2 & 
1.5 for Lya LFs at z —3.1, 4.5, 5.7, and 6.6, respectively. 

5. RESULTS 

Our model Lya LFs with single star formation effi- 
ciency (SEE), agree remarkably well w ith the observa- 
tions. To predict the Lya LFs of LAEs. lNagamine et al.l 
(|2008l ) investigated two models, the duty cycle & escape 
fraction scenario, and found that the duty cycle model 
reproduces observations better than the escape fraction 
model. In our model, the duty cycle is naturally pro- 
duced since only halos with high accretion rates will be 
observed as LAEs. Fig. 2 shows the halo mass function at 
z = 4.5 (blue dashed line) and the accreted mass function 
(blue solid line). Thus, the use of accreted mass rather 
than the total halo mass eliminates the need to introduce 
an additional duty cycle parameter in our model. 

A constant SEE over a wide redshift range implies that 
the physical properties including Lya luminosity, SFRs, 
and stellar masses of LAEs do not evolve significantly 
from z » 3 — 7 sinc e these properties depend on SEE. 
iJimenez et al.l ()2005l ) also found a similar result of con- 
stant SEE (converting accreted baryons to stars) over a 
wide range (w 2 orders of magnitude) of stellar masses, 
and over a relatively large redshift range, using a large 
SDSS sample combined with stellar population models. 
Fig. 3 shows the DM accretion rate as a function of halo 
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mass. The solid line is the least-square fit to the median 
mass in a AM/dt bin. A nearly constant (0.8-0.9) slope 
of this line at all redshifts implies that the SFR does 
not evolve in this redshift range, while deviation of the 
slopes from unity suggests that the mass accretion rate 
is a non-linear function of halo mass. 

5.1. Physical properties of LAEs 

We use the best-fit model parameter (SFE) at z = 3.1 
to derive other physical properties of LAEs and compare 
our results with the observations. Our best-fit model 
yields an SFE of 2. 5%, c onsistent w ith th e global star for- 
mation effici ency. iFuku gita ct al| (119981 ) predict a SFE 
< 5%, while iBaldrv et al. ( 2008.) found this value in the 
range 4 — 8% for blue light in galaxies. 

The SFE of 2.5% yields SFRs « 1 - 10 Mq yr"! 
corresponding to the observed Lya luminosity range 
LLya ~ 1 X 10*2 _ I X 10^^ gj.g ^-1 ^ ^his SFR is Com- 
parable to the inferred SFR « 8 Mq yr^^ in LAEs at 
z w 5 (jPirzkal et al.ll2007t ). A similar average value 
of SFR pa 6 M(:) yr ~^ was inferred for z — 3.1 LAEs 
(|Gawiser et al.ll2006[) . A sUghtly higher value of SFR 
« 5.7 - 28.3 Afgyr-i with med ian SFR ft! 9.6 M(^yr~ ^ 
was inferred for z = 5.7 LAEs ([ Muraya ma et al.l 120071 ). 
For z = 6.6 LAEs, iTaniguchi et al. (200^ found an av- 
erage SFR « 5.7 ±2.3 M0yr~^. These averages however 
depend on the depth of the surveys; deeper surveys probe 
less luminous galaxies and hence lower SFRs. 

The total stellar mass in young stars (estimated us- 
ing equation 3) of LAEs corresponding to the ob- 
served Lya luminosity range is M^, « 3 x 10^ — 3 x 
IO^Mq, in good agreement with the observed stel- 
lar masses of lO'^'Mf;^ — lO^Mp, for a LAE sample at 
z = 4.5 fFinkelstein et all l2007t iGawiser et al] l2007l : 
iPirzkal et al.l 120071 : iPentericci et al.l l2009f ). Thus, our 
model reproduces the primary physical properties of 
LAEs at z « 3 — 7. In addition, our assumption of large 
escape fraction of Lya photons will yield high observed 
equivalent widths of Lya line in LAEs. 

5.2. Evolution of Lya luminosity function 

The Lya LFs have been used to probe the epoch 
of reionization and constrain the evolution of in- 
tergalactic medium (IGM) (e.g. Haiman & SDaans" 
1999; Malhotra & Rhoads 2004; Haiman & Gen 2005; 



Stern et al.l 20051: iK ashikawa ct al. 2006; Diikstra ct al 
2007HMcQuinn et al..,2007; Mesinger fc FurlanettQ,200S : 



Ota et al. I2008D . Any significant evolution of the Lya 



LF, after correcting for the intrinsic change in the num- 
ber density of LAEs, between two redshifts will imply 
that the IGM evolved at these redshifts. 

Currently the evolution of Lya LF at z > 5 is not 
well understood. Previous studies find no significant 
evolution in Lya LFs at redshifts between z = 5.7 and 
z = 6.5 (jMalhotra & Rhoads 20 l ). However rece nt ob- 
servations (IKashikawa et al.ll2006rfota et al.ll2008D find a 
modest decline in the bright end of the Lya LF from z = 
5.7 to z = 6.6, suggesting IGM evolution at these red- 
shifts. iDiikstra etall (|2007D showed that the weak IGM 
evolution between z = 5.7 & z = 6.6 can be attributed 
to the evolution of DM halo mass function. In addition, 
the cosmic variance in a volume limi ted LAE survey also 
affects the Lya LFs. For example, IShiova et al.l ()2009D 
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Fig. 4. — Field-to-field variance of number of LAEs at 2 = 6.6, 
measured in eight subvolumes (102x51x25 Mpc^), plotted as 
a function of survey detection limit. The variance > 30% 
for a typical narrow-band LAE survey with Lyo detection limit 
> 2 X 10"^ erg s'^. 

find that the number density of LAEs, at z » 5 vary by 
a factor » 2 in a survey area of 0.5° x 0.5°. We now 
investigate the evolution of Lya LFs, and the effect of 
cosmic variance on number density of LAEs in a volume 
and flux limited LAE survey. 

From Fig. 2 we see that there is some evolution of num- 
ber density of DM halos that can host LAEs. However, 
this weak evolution is due to the intrinsic change in the 
number density since we haven't included any IGM cor- 
rection in our model. Thus our model can explain the 
observed evolution of Lya LF (jKashikawa et al.l l2006l ) 
without invoking the cosmic variance. However, it would 
be interesting to estimate the effect of cosmic variance on 
the Lya LF. 

In order to understand the effect of cosmic variance 
on the observed number density of LAEs in a volume 
and flux limited LAE survey, we estimate the field-to- 
field variance by dividing the total simulation volume 
into eight non-overlapping rectangular boxes, each with 
a comoving volume (102x51x25 Mpc^) comparable to 
typical narrow-band LAE surveys. The variance is cal- 
culated using —< {N — /i)^ > //^, where N is the 
number of LAEs in each sub- volume and /i is the aver- 
age number of LAEs. Fig. 4 shows field-to-field variance, 
for a volume limited LAE survey with different Lya flux 
limit. For a narrow-band LAE survey with Lya detec- 
tion Hmit oiL^ya > 2 x 10^^ erg s~^ and a survey volume 
< 2 X 10^ Mpc'^, the field-to-field variation is significant 
with > 30%. This result confirms the necessity of us- 
ing a large volume to minimize cosmic variance in LAE 
surveys. It also strengthens our conclusion that some 
apparent evolution of Lya LF can be attributed to the 
cosmic variance. 

6. CLUSTERING OF LAES 

LAE s are found to t race rarer, highe r density re- 
gions (|Iliev et all 120081 : lOrsi et al.l l2008f ). However, 
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Fig. 5. — The spatial distribution of our model LAEs in a slice from DM simulation at two redshifts z = 5.7 (left) and z = 4.5 (right) 
in a volume 30 X 30 X 17 Mpc^. The small (red) and big (blue) filed circles represent the positions of DM halos and model LAEs, 
respectively. Only LAEs with Lj^ya > 2 X 10*^ erg s~^ are plotted. In general, the LAEs are located in high density regions. Also note that 
different halos host LAEs at the two redshifts, depending on whether they are accreting or not. This gives rise to a duty cycle quite naturally. 



IShimizu efall (|2007t ) suggested that the LAEs at z « 3 
do not necessarily reside in high density peaks. In this 
section, we investigate whether our model LAEs reside 
in high density regions, and estimate their spatial corre- 
lation lengths. Fig. 5 shows the spatial distribution of 
LAEs at two different redshifts, z = 5.7 and z = 4.5, in 
a simulation slice of 30x30x17 h^^ Mpc^. The depth of 
this slice is comparable to the depth of a typical LAE sur- 
vey. Comparing the locations of LAEs at two redshifts 
it is clear that different halos host LA Es at different red- 
shifts, thus exhibiting a duty cycle (jStark et al.l l2007t 
iNagamine et al.ll2008l) . Our model LAEs are generally 
located around overdense regions consistent with the ob- 
servations and as expected in biased galaxy formation 
models . 

6.L Two-point spatial correlation function 

The two-po int spatial correlation function i^(r) 
(jPeeblesI Il980f ) is frequently used to study the clus- 
tering properties o f galaxies (e.g. IZehavi et all 120051 : 
iGawiser et al] 120071: iKovac et al. l l2007f). We use the 
Landy -Szalay estimator, proposed bv iLandv fc Szalavl 
(|1993f ). to calculate the two-point spatial correlation 
function given by 

CM 



DD{r) - 2DR{r) + RR{r) 



(5) 



RR{r) 

where DD{r), RR{r), and DR{r) are the number 
of galaxy-galaxy, random-random and galaxy-random 
pairs, respectively with separation distance of {r,r + 6r). 
To compare our model ^(r) with the observations and 
quantify its evolution with redshift we only include LAEs 
brighter than the detection limit of LAE surveys at 
a given redshift, and use our full simulation volume 
w 1 X 10^ Mpc'^ for better statistical significance. 



To calculate ^(r) at each redshift, we generated a ran- 
dom sample of points with uniform coordinates drawn 
from a uniform probability distribution, and a number 
of random points exactly equal to the number of LAEs. 
We count the number of pairs, DD(r), RR(r), and DR(r) 
separated by a distance r by binning the points at differ- 
ent r with binwidth of 5r = 0.2h~^ Mpc. To minimize 
the random errors, we perform 50 realizations with dif- 
ferent sets of random points and calculate an average 
^(r) at each r. Using r < 20 Mpc and assuming 
negligible error, we obtain the spatial correlation length 
tq by fitting a least-square power law to the correlation 
function. 

The correlation lengths obtained from our model LAEs 
show a modest evolution, with ro = (3.2 ± 0.3, 5.0 ± 
0.3,4.2±0.6,6.0± at z =(3.1, 4.5, 5.7, 6.6), re- 

spectively. Fig. 6 {left) shows a comparison between our 
model predictions (filled circles) and observed tq (shown 
with different symbols given in the labels) at different 
redshifts. The model predicted ro values are slightly 
shifted along x-axis to avoid overlap with the observa- 
tions. The predicted correlation lengths are consistent 
with the observations at z = 3.1 with observed ro = 2.6 
±1 (jGawiser et al.ll2007D. and at z = 4 .5 with ob- 
served ro = 4.6 ±0.6 (|Kovac et al.|[2007l ) estimated 
using contamination-corrected (the maxim um value per- 
mitted) LAE sample. However. I Ouchi et a l. (2003) found 
a higher ro = 6.2 ± 0.5 Mpc for contamination- 
corrected LAE sample z = 4.86. It would be interest- 
ing to estimate the variance of ro in a volume and flux 
limited LAE survey, and see if we can account for the 
large difference seen in the observed ro between z =4.5 
& z =4.86 LAE surveys. 

To estimate the variance of ro (c^o) ^ volume and 
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Fig. 6. — Correlation lengths of LAEs at different redshifts. Left. Comparison of correlation lengths of our model and observed LAEs 
at different redshifts. Here we include LAEs with Lya luminosity greater than the survey limit at each redshift. The filled circles are our 
model results, other symbols are from different observations as shown in the labels. The observed ro values shown here for z=4.5 & 4.86 
are for the contamination corrected (the maximum value permitted) LAE sample. Our model results are slightly shifted to avoid overlap 
with other observational points. Right. The correlation lengths of our model LAEs with a constant Lya luminosity cutoff (Li^y^ > 1 X 10^^ 
erg s~^) at all redshifts. 



flux limited LAE survey, we divide the total simulation 
volume at z =3.1 into five non-overlapping sub-volumes, 
each with a comoving volume of (102 x 102 x 20.4 Mpc'^), 
approximately equal to a typical survey volume and only 
including LAEs with L^y^ > 1 x 10^^ erg s~^. We calcu- 
late ^(r) and rg, as described above (second paragraph) 
in each sub- volume and estimate a"^^^ at z =3.1. We find 
(7^^=0.5 Mpc with an average ro=3.2 Mpc, av- 
erage of To in five sub- volumes. While this variance in 
Tg cannot account for the large difference in rg observed 
between the two surveys at z =4.5 & 4.86, it is clear 
that one needs to take into account such variance in cor- 
relation lengths obtained from volume and flux limited 
surveys. 

We now investigate if the redshift evolution of rg seen 
in Fig. 6 {left panel) is significant since this redshift 
evolution can result from the surveys at lower redshifts 
extending to lower luminosities and hence probing lower 
halo masses. In order to understand this effect, we con- 
sider full simulation volume and only include LAEs with 
a constant L^ya > 1 x lO**^ erg s~^ at all redshifts. 
Choosing a constant lower luminosity at all redshifts 
implies that we are probing approximately same halo 
masses at all redshifts. We calculate rg in the same way 
as described above (second paragraph) except that we 
impose a constant luminosity cutoff at all redshifts. Fig. 
6 {right panel) shows rg with a constant lower luminosity 
at all redshifts. While the predicted correlation length 
show some evolution with redshift, even after applying a 
constant lower luminosity cutoff at all redshifts, most of 
this evolution seen in Fig. 6 {left panel) can be attributed 



to the luminosity limit of different surveys probing differ- 
ent halo masses. In order to test our model predictions, 
we suggest that more LAE observations are necessary for 
probing a wider range of halo masses especially at z >5. 

7. SUMMARY AND CONCLUSIONS 

We have used a physical model with a single vari- 
able parameter to populate dark matter (DM) halos 
with LAEs in a cosmological simulation and compared 
our model predictions with the observations at redshifts 
z « 3 — 7. In our model we assumed that the star for- 
mation rate (SFR), and hence the Lya line luminosity 
is proportional to the mass accreted by halos. In other 
words, the star formation in LAEs mainly results from 
the cold accretion of new material. Despite the lack of 
observational evidence relating accretion rate to the Lya 
luminosity of LAEs, it is promising that our model gives 
good fit to the observations over a wide range of redshifts 
and is able to reproduce several other physical proper- 
ties of LAEs. In addition, relating the accreted mass 
(rather than the total halo mass) to the Lya luminosity , 
gives rise to a duty cycle of LAEs (jNagamine et al.ll2008f ) 
quite naturally. 

To compare our model predictions with the observa- 
tions, we first constructed the Lya luminosity function 
(LF) at z = 3.1 and obtained the best-fit model by vary- 
ing the star formation efficiency (SFE) and comparing 
the model Lya LF with the observations at this red- 
shift. We then used this best-fit model to predict the 
Lya LFs and physical properties of LAEs at z = 4.5, 5.7 
and z — 6.6. 

Using a constant SFE, our model predicted Lya LFs 
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agree remarkably well with the LAE observations over 
a wide redshift range. Since other physical properties 
including SFR and stellar masses of LAEs depend on 
SEE, it implies that these properties do not evolve sig- 
nificantly over these redshift range. Our best-fit model 
yields a SEE=2.5% which gives SFR « 1 - 10 MQyr~^ 
in good agreement with the observations. We find that 
the model LAEs in the currently observable luminosity 
range (2 x 10"*^ < L^ya < 2 x lO'^'^erg s~^) have stellar 
masses sa 3x 10^ to 3x lO^M© of young (< 30 Myr) stars. 
These stellar masse s of LAEs are similar to those inferred 
from observations (|Finkelstein et alll2007t iPirzkal et al.l 
|2007( ). We also investigated the evolution of Lya LEs 
from z « 3 — 7 and find that there is no significant evolu- 
tion of Lya LE due to the IGM if we include the intrinsic 
change in the number density of LAEs over this redshift 
range. This conclusion is strengthened if we include the 
effect of cosmic variance on the observed number density 
of LAEs. We show that the field-to-field variance can 
be large ~ 30% for a flux and volume limited surveys 
comparable to current observations. 

We studied the clustering properties of LAEs and 
found that the LAEs are mostly located in the high 
density peaks. Our model predicted correlation lengths 
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